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We have studied the dissipative dynamics of a driven electronic spin trapped in a quantum dot. 
We consider the dissipative mechanism as due to the indirect coupling of the electronic spin to 
acoustic phonons via the spin-orbit /electron-phonon couplings. Using an effective spectral function 
of the dissipative phonon bath, we evaluated the expectation values of the spin components through 
the Bloch-Redfield theory. We show that due to a sharp bath resonance present in the effective 
spectral function, with typical energy much smaller than the electronic confinement energy, the 
dissipative spin has a rich dynamical behavior that helps us to determine some features of the spin- 
bath coupling. We also quantify the effects produced by the sharp bath resonance, and thus indicate 
the best regimes of operation in order to achieve the longest relaxation times for the spin. 



I. INTRODUCTION 

During the last decade several physical systems have 
been proposed as candidates for qubits- and, among 
those, condensed matter devices have been attracting a 
growing interest in the area 2 -. There are many reasons for 
this and most of them involve the way one designs, con- 
trols and accesses individual qubits in quantum systems 
containing a very large number of such entities. This is 
a natural consequence of the fact that dealing with solid 
state devices one can use the same standard electronic 
circuitry of conventional computers. Furthermore, new 
and powerful experimental techniques allow us to expect 
that these qubits could be more easily scalable and that 
control over the processing of information and implemen- 
tation of possible protocols of error correction could be 
done in a more reliable and efficient way. These are only 
a few points in favour of these systems as being good 
candidates for qubits. 

Despite these important favourable points, there are 
also negative aspects concerning the use of solid state de- 
vices as qubits. Since they are mainly meso or nanoscopic 
devices, which under specific conditions mimic two-state 
system behavior, it is a very hard task to isolate them 
from their environment. Moreover, as it is already 
well known, the system-bath interaction causes loss of 
quantum coherence which is a drawback for quantum 
computation^, in particular, if we do not wish to op- 
erate our devices at too low temperatures. Therefore, 
there must be a compromise between the desirable fea- 
tures possessed by these systems and the undesirable ef- 
fect of decoherence. It is along this direction that people 
have been investigating electronic spins in quantum dots 
as promising candidates for qubits^. Since the confine- 
ment of the electronic wave function isolates the spins 
from most energy relaxation channels, there are situa- 
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tions where decoherence takes place only beyond a fairly 
long time interval. Here we mean a time interval within 
which a large number of logical operations could be coher- 
ently performed. In particular, for self-assembled quan- 
tum dots, the energies involved in the spin dynamics are 
such that the decoherence time is indeed very long. How- 
ever, what is advantageous on the one hand turns out to 
be a problem otherwise. As the spin is a microscopic 
variable one has to face the problem of how to access it 
experimentally. 

A possible way to address one of the many spins in 
the quantum processor is by tuning the frequency of a 
time dependent external magnetic field (the transverse 
pumping or control field) to the spin's Larmor frequency. 
Indeed, Koppens and collaborators 7 have demonstrated 
very recently the first experimental realization of single 
electron spin rotations in quantum dot systems, which 
is a necessary step for the implementation of universal 
quantum operations. They showed, through current mea- 
surements in double quantum dots, the ability of per- 
forming spin-flips using a resonant oscillating magnetic 
field with the spin's Larmor frequency. Therefore, in this 
work, we shall investigate the dissipative dynamics of an 
electronic spin confined to a quantum dot subject to a 
strong static magnetic field and a much weaker trans- 
verse pumping field. 

It has been demonstrated that the main mechanisms 
of relaxation of the latter are the spin-orbit interaction^ 
and the hyperfine interaction with the spins of the host 
lattice^. For static fields stronger than IT it is expected 
that the spin-orbit interaction plays a major role in the 
relaxation process. The relaxation times experimentally 
observed for fields above 4T are within the range ~ 0.8 — 
20m a 11 i 12 which encourages us to study the possibility of 
implementing spins confined in quantum dots as qubits. 

Following previous worksi^, we investigate the dissipa- 
tive effects originating from the coupling of the orbital 
electronic motion to the acoustic phonon modes of the 
lattice where the dots are embedded. Due to the spin- 
orbit coupling, the spin degree of freedom becomes indi- 
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rectly damped by the latter. Consequently, it is possi- 
ble to define an effective dissipative two-state system^ 
whose dynamics may not be describable pertubatively-^. 
In particular, one should stress the existence of a very 
sharp effective bath resonance at energies much lower 
than the planar electronic confinement energy. 

The main goal of the present work is to evaluate the 
dissipative dynamics of a structured spin-boson model in 
the presence of a time dependent magnetic field aiming 
at the determination of the best operational conditions 
under which these systems could be employed as good 
qubits. In section [TT] we present the model we use to 
describe this dynamical process. In section 1111} we ob- 
tain the approximate analytical solutions to the Bloch- 
Redfield equations for the average values of the spin com- 
ponents. Finally, in section Hvl we carry out a detailed 
analysis of the solutions we have obtained. 
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Figure 1: Spectral density J e // for piezoelectric electron- 
phonon interaction, s — 3. Inset: The region close to the 
peak fi 3 can be well approximated by the Lorentzian in Eq. 
[5] (dashed line). Here we assumed ■^ a - » 0.5 and 5$ — 355. 



II. MODEL 

We consider quantum dots with strong confinement in 
the z direction, and electronic orbital motion in the x — y 
plane subject to a confining parabolic potential. The 
spin-orbit coupling is modelled by a Dresselhaus interac- 
tion term projected onto the plane of the dot. Besides, 
we add a term due the presence of an externally applied 
magnetic field B. Thus, these assumptions lead to the 
following spin-orbit Hamiltonian 



H. 



so 



hiuo (a\.a x + - l3a x P x 



+hu) [a\a y + ^ + PvyPy + t^bB • <r, (1) 

where ujo is the lateral harmonic frequency, f3 = j c (k 2 ) = 
r y c m*uj± (with j c is the Kane parameter—, to* is the elec- 
tron effective mass), or are the Pauli matrices, and a x (y) 
is the usual ladder operators for the x(y) direction. 

The electron-phonon Hamiltonian, considering either 
the piezoelectric or the deformation potential interac- 
tions with acoustic phonon modes, can be mapped into 
the bath of oscillators model?! with the spectral function 
given by! 3 - 



J s (w) 



m*u) 2 D 8 s 



CO 
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where s — 3 for the piezoelectric interaction with dimen- 
sionless coupling S 3 = { ^$f + js) , and s = 5 



for the deformation potential with 5$ = 



a cr"c w here 

cud is the Debye frequency, vi and vt are the longitudi- 
nal and transverse sound velocities respectively, p is the 
material density, (e m )i4 is the electromechanical tensor 
for zinc-blende structures^, and a Ci r is the deformation 
potential at the T point^. 9 is the Heaviside step func- 
tion. 



Provided that we have introduced a bath of oscilla- 
tors for the electron-phonon coupling, and because of the 
spin-orbit interaction, one can think of this problem as 
a spin degree of freedom coupled to an effective bath 
of oscillators^ with the following effective bath spectral 
function 



j eff H = TO*/? 2 



s+2 



-6 {uj d - to) , (3) 



Z{u) +S* l£ 



where Z( W ) ee - [^f (l + S s c/> s (^)) , and 

0.0c) = fP/oV/O/ 3 -yx 2 )dy = -2^1(B(x,8,0) + 
(— l) s B(— x, s, 0)), with B being the generalized incom- 
plete beta function. 

Figure Q] presents the J e //(w) behavior for the piezo- 
electric interaction case, s = 3. The most prominent 
feature of its behavior is the peak at u> — fl s , where f2 s 
is determined by 



(,s-2) 
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Furthermore, the region around fl s can be 
well approximated by a Lorentzian of width 

s-1 



ure[T]), given by 
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(see inset of Fig- 
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The asymptotic analysis of J e //(w) also reveals that in 
the low frequency range, defined by uo <C il s and <C 

( u>t T ) > tne effective spectral density is always super- 
ohmic, with power s + 2. On the opposite side of the 
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resonance, and in the high frequency limit, l! s « w < 
iv d, the spectral density behaves as a power law in the 
frequency w with exponent s — 2. 

Since we know the spectral density of the effective bath 
of oscillators coupled to the spin degree of freedom, our 
problem can now be modelled by the driven spin-boson 
Hamiltonian 



-\Aa x + e(t) a z 



where the spectral density of the environment is 
given by |J3J, and the applied magnetic field B = 
(— HA/gnB,0, — tve{t)/gfiB)- In the weak coupling limit, 
the equations of motion for the spin expectation val- 
ues, <Tj (t) = Tr Ip (t) can be written, for this specific 
model, as the generalized Bloch-Redfield equation 16 

& x (t) = e (t) <j y - T xx (t) a x - T xz (t) a z - A x (t) , 

&V = ~ £ °X + A(7 2 - Tyy (t) (Jy ~ Ty Z (t) (T Z ~ Ay (t) , 

& z (t)= -Aa v , (7) 

where the fluctuating terms are given by^ 

A x (t) = ImF (t) , A y (t) = ReF (t) , with 

F(t) = 2 f dt'lmM(t - t')U RR (t,t')U RL (t,t'), 
Jo 

and the temperature dependent relaxation rates are de- 
termined by-i£ 



r I3 (t) 



dt'Re[M(t - t')]bij(t,t'), with F yy (t) = Y xx (t) . 



The correlation function 
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M (t) = - duj.J(uj) 



cosh(^ - iujt) 



sinh(^) 



with j3 = 1/ksT, takes into account all the effects of the 
bath of oscillators. The functions fry read 

b xz = 2Re [Ur R U rl ] , b yz = -21m [U RR U RL ] 

and b xx = \U RR \ 2 - \U RL \ 2 , where U(t,t') is the non- 
dissipative time evolution operator for the spin system 
and its matrix elements U RR (t,t') = (R\U(t,t')\R) and 
U R L(t 7 t') = (R\U(t,t')\L), are written in the basis de- 
fined by a z \R) = + \R),a z \L) = -\L). Eqs. © are 
derived^ under the assumption that : ^f L In(f2 c /A) <C 1, 
where Sl c represents the cutoff frequency of the system. 



considering it as a small perturbation to the spin dynam- 
ics, i.e., e /A <C 1. In order to solve the Bloch-Redfield 
equations we need to compute the non-dissipative time 
evolution operator U(t,t'). Using the rotating wave ap- 
proximation (RWA), we obtain the following simple an- 
alytic form to the time evolution operator 



U(t,t ) = Rl (fii)e 



i€l(t—t )c 



•ijxCmo), (8) 



e 2 



A-il^ + ^a z 



where R x (fit) 

cos ((f)) a z + sin (cf>) a x , with tan (0) = A ~ n and t\ = 

e§ + (A — SI) 2 , ei and A are, respectively, the Rabi 
and Larmor frequencies of the problem, SI is the exter- 
nal field frequency and 4> represents the frequency de- 
tuning. The Bloch-Redfield coefficients can now be eval- 
uated using the spectral density of the bath oscillators 
and U (t,t r ). We proceed evaluating the time integrals 
first. Those have sine and co-sine elementary integrands 
that, when appropriately rearranged, can be written in 
terms of sine functions, whose maxima occur at the nat- 
ural frequencies of the system SI, and SI ± ei. Follow- 
ing this procedure, the coefficients of the Bloch-Redfield 
equations can be written in the compact form 



A x (t) 

Ay (t) 

r xx (t) 
r xz (t) 
r„* (t) 



n (t) + mi (t) sin (2Slt) + m (t) cos (2Slt) , 

77i2 (t) sin (fit) + Tl2 (t) COS (Sit) , 

r 3 (t) + TO 3 (*) sin (2Slt) + n 3 (t) cos (2Slt) (9) 

7774 (t) Sin (fit) + 774 (t) COS (Sit) , 

r 5 (t) + 7775 (t) sin (2Slt) + n 5 (t) cos (2Slt) , 



(see appendix for the explicit presentation of those co- 
efficients). The new coefficients Ti (t) , m, (t) and m (t) 
have their time dependence determined by the integrals 
I n , Eqs. (fT9ll and (J20J) - Those integrals, in the regime 
of interest, t 3> l/u>u, have their main contributions 
decomposed in two parts: one, time-dependent, occurs 
due to terms arising from the poles of the spectral den- 
sity in the complex plane, which have lifetimes deter- 
mined by their imaginary parts; the other part is a time- 
independent contribution due to the resonance of system 
with the external applied field. Therefore, each coeffi- 
cient Ti (t) , mi (t) and 77^ (t) can also be written in well- 
characterized time dependent and independent parts. In 
fact, considering J(uj) given by ([3]), whose main poles are 
fl s ± it, we can write n (t) as 



n(t) ^U^e-^M V-flW 



i=0 k 
sin[(O s + fi (i) )t] 
Sl s + Sl^) 



■n, 



(10) 



III. RESULTS 

From here onwards we will be interested in the study 
of a monochromatic field of the form e (t) = 2e<j cos Sit, 



where Sl^ = SI, Sl^ = SI + ei and Sl^ = SI — ei; 



Co = cosW/2, Ci 



(l + sin(0))V4, and C 2 



(1 - sin(0)) /4. The first term on the r.h.s of Eq. (10} 
is a direct contribution of the poles of J (u>) , and the 
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last term, fi, takes into account the time independent 
terms of r\ (t). Substituting ifTO)) in the Bloch-Redfield 
equations, we can find a particular solution for a x , up to 
first order in £, due to the first term of Eq. (flO| . Thus 
we can write a x (t) as follows 



a x {t) a £J(Q s )e- it J2C i 



(1 - cos (n s - nM) t) 



l=0 k (n, - n«>) 

(1 - cos + fiW) i) ' 



5-x (i) • 



It is clear that the particular solution has a lifetime deter- 
mined by the width I of the peak of the effective spectral 
density. The modes arising from this term oscillate with 
the system's natural frequencies shifted by the bath res- 
onance f2 s . Thus, the contributions due to the poles of 
J e ff(ui) introduce a new time scale, 1/1, in the dynamics 



of the dissipative spin, which distinguishes the short and 
long time regimes of the problem. The first term of lflO|) 
is expected to be the dominant contribution of the poles 
of J e ff(uj) for the spin dynamics, because those arising 
from mi(t) and n\(t) oscillate with frequency 20, larger 
than those from n (t) . After integrating the equations of 
motion over any measurable time interval, their effects 
become negligible. 

Now, since we have determined the coefficients 
r(t), m(t), and n(t), we can obtain analytic solutions for 
the Bloch-Redfield equations. For the long time regime, 
t > l/£, those coefficients approach very quickly their 
asymptotic values f = f(t —* oo), with / = r, m, and n. 
A solution can be obtained using the Laplace transfor- 
mation for the regime of interest, namely eo/A <C 1. 
Retaining terms up to order of £q, and considering the 
initial condition <r x (0) = ±l,a y (0) — 0,er 2 (0) = 0, we 
find the following approximated solutions 



<J X (t) 



U (fi.) e- H ]T a I V X om ; ' + V X - om 1 1 } + a x (oo) + G cos (2fit) 



i=0 



2 \ 



+ o x (0) - a x (oo) - G - Hi ) e- r « ( + £ i^e"^* cos (<**) , 

\ i=0 / i=0 

^/ (e Q , 0) (1 - cos (Qt)) | ^ (0) -a x (oo)-G-J2 e ~ rRt + (°°) 
+A £ { f ^ n -^ (i _ C os ((17 - Ui ) t)) + + (1 - cos (( i I ) / i ) 



r 



(ii) 



(12) 



where 



LU2 



g+(fl- v/A(A-2 7 )) « ^ + (tt-A + 7 )) 2 



eg + fn+^/A (A-2 7 ) «*/«=§ + (n + a - 7 ) 



The shift 7 in the natural resonance frequency can be un- 
derstood as the frequency correction to the spin dynamics 
due to the interaction with the bath of oscillators. In con- 
trast to the procedure adopted inAS, this correction arises 
naturally from the solution obtained. The frequency shift 
7 has the form of a Lamb shift 



7 ~ r 5 



E 

i=0 



Ci 
2t Jo 



Wd ( Huj 

duj coth ( ) J (w) P 



\2k B T 



(13) 



r 



This expression is very similar to that found in the non- 
perturbative treatment of the atom-field interaction^. 
However, it is worth noting that Eq. (fl3|) is derived here 
as a temperature dependent function. 

The rates appearing in Eq. ifTTj) and (fl2l) can be writ- 



ten approximately as T R « ^1 - ^/^"'(Tn) A( 'n^ 27 ' > ' ' 
and Tq = 2Ti = = 2T , where 



r ~3 = Iyc ! j(oW)cothf-^ 

z=0 v 



(14) 
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Figure 2: Rate F, Eq. I14I as a function of the driven field 
frequency fl (solid line). There is a very pronounced peak at 
the region close to f2 s . Dashed line represents the rate F for 
the non-driven case, eo = 0. Temperature used T fs IK. 



Figure [2] presents the rate T as a function of the driven 
field frequency fi. It is clear that there is a signature 
of the peak of spectral density J e // (w) in its behavior. 
Close to the peak, O ps fi s , the rate T can be several 
orders of magnitude higher than the asymptotic value 
\ J (A) coth (fi,A/2fcsT). This result shows which regime 
of frequency must be avoided to obtain the largest relax- 
ation time of the system. 

Finally, the coefficients can be written as 



<j x (oo) 
G 

Ho 



l + 2/ 2 (e ,fi) 



tanh 



/ HA 



^/ (e ,O)o- x (oo) , 



-2G, #i 



<J X (0) . 



The function 



/(eo,fi) = eo 



fi 



fi 2 



A (A - 2 7 ) ' 



gives a measure of the effects of the driven field on the 
spin dynamics. As expected, there are two important 
features to be considered: one is its intensity heo/g/iB', 
and the other is how di stant fi is fr om the shifted natural 
resonance frequency -y/A (A — 27) « A — 7. 

In the limit eo 0, we obtain a x (t) ps a x (00) + 
(cs (0) — <J X (00)) e~ rt , where T is the well-known 



expression 
tanh' hA 



T = iJ(A)coth(fiA/2fc B T), and 0^(00) = 
is exactly the expected thermodynami- 

For the resonant case ft ~ A — 7, 
- 7) — > 00, which implies a x (t) sa 
The system evolves with a frequency 



2k B T / 

cal mean value, 
we have /(eo, A - 
o- x (0)e~ rt cos(e t). 
imposed by the external driving field, and relaxes to a 
state where both spin components are equally probable. 
A detailed description of those phenomena is presented 
in the next section. 
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Figure 3: Time evolution of the expectation value of a x for the 
near resonance case Q ~ A — 7. First plot presents the initial 
dynamics, where two well defined oscillation frequencies can 
be seen. For times t < l/£, damped oscillations due the poles 
of J e ff(uj), with relaxation time l/£, are observed. In the 
long time regime, t > l/£, the system evolves with oscillation 
frequency eo and, initially, with large amplitudes (reflecting 
the tuning between the frequency of the driving field and the 
natural frequency of the system). Due to the dissipative pro- 
cess, the system tends, in the time scale 1/T, to an equal 
probability state for spin up and down, vanishing the average 
of a x (second plot). For very long times, t ^> the system 
reaches a stationary state, oscillating with frequency 2f2 and 
very small amplitude. Here we assumed loo/£ = 750. Solid 
and dashed curves represent the numeric and the analytic re- 
sults obtained. 



IV. DISCUSSION AND CONCLUSIONS 

Figures [MS] present the expectation values of a x and a z 
for the near resonance, ft « A — 7, and the off-resonance, 
Cl = A, cases. Here we assume hA ^ fceT, A « fig/2, 
e /A fa 0.03 and uj d /£ ps 750. 

It is well knowni^ that, in the absence of dissipation, 
the system would indefinitely evolve into coherent cycles 
of absorption and emission of energy - the Rabi oscil- 
lations. Furthermore, the change in the spin up and 
down populations exhib its an oscillato ry time depen- 
dence with frequency \J (\ + (fi — A) 2 , and amplitude 
e o/( e o + (fi — A) 2 )- It is worth noting that under the res- 
onance condition, fi = A, the oscillating field fully drives 
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Figure 4: Time evolution of the expectation value of a x for 
the off-resonance Q. = A case . Likewise the resonant case, 
the initial dynamics is composed of two well denned oscil- 
lation frequencies (first plot). Again, the contributions due 
the poles of J e ff(u) vanish after times t « If I. In the long 
time regime, t > l/£, the system evolves oscillating with fre- 
quency « \J £q + 7 2 , and the amplitude observed is not as 
larger as in the resonant case. The second plot presents the 
dynamics for the long time regime. As one can see, the sys- 
tem tends to approach the stationary state at the same time 
scale as the resonance case, however we see that a x (t — > oo) 
is now much closer to its thermodynamical equilibrium value 
for eo = 0, t&nh(hA /2ksT) ^ 1. Once again, for very long 
times, t 3> 1/F, (inset second plot) the system oscillates with 
frequency 2Q, and very small amplitude. Here we assumed 
lod/£ = 750. Solid and dashed curves represent the numeric 
and the analytic results obtained. 



the transitions between up and down spin states and that 
the weaker the driving field, the narrower the resonance 
peak. As one can notice by inspecting figures [3] and HI 
the coupling system-bath changes some of those features. 
First of all, the resonance condition is no longer verified 
at the frequency associated with the static applied field. 
Now, the new natural frequency is given by fl = A — 7, 
where the shift 7, Eq. lfl3|) . is completely associated with 
the coupling between the system and its reservoir - Lamb 
shift. This is exactly the behavior we have for the first 
plots of figures [3] and HJ the largest amplitudes of os- 
cillation occur for those frequencies closest to the new 
resonance condition. As expected, the component a x os- 
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Figure 5: Time evolution of the expectation value of a z for 
the near resonance case SI ~ A — 7. The initial a z dynam- 
ics behavior (first plot) is characterized by a clear structure 
of beats, with beat and angular frequency given by uii and 
Q, + LUi/2. Due the decoherence processes, the modes related 
with the frequencies ua evanesce in the characteristic time 
1/r, destroying the beat structure (second plot), for which 
arises an oscillatory regime of same frequency f2 of the applied 
field. Here we assumed u)r>/£ = 750. Solid and dashed curves 
represent the numeric and the analytic results obtained. 



dilates with frequency eo at resonance, and \Ao + 7 2 f° r 
the off-resonance fl — A case. Another important change 
in the spin dynamics is that the system does not evolve 
indefinitely in coherent cycles of absorption and emission 
of energy. In fact, the dissipation process induced by 
the coupling to the dissipative bath destroys the coher- 
ence in the emission-absorption cycles and, after a long 
time, decoherence takes place. Indeed, the behavior of 
the expectation value of a x (t) in the long time regime 
(figures[3]and[4|), reveals this decoherence process. As one 
can see, the large amplitude of oscillations observed ini- 
tially, a consequence of the coherent emission-absorption 
stimulated by the field, tends to decrease in time. This 
process occurs in a time scale given by 1/r, Eq. ffl4|l . 
which has almost the same value as presented in figures 
[3] and [H This happens because for both cases we are 
away from the peak of the rate T, Figure El where it is 
practically a constant given by \ J(A) coth(hA/2kBT). 
For very long times, t 1/r, the system reaches a sta- 
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Figure 6: Time evolution of the expectation value of a z for 
the off-resonance case f2 « A. The same behavior for the res- 
onant case is observed. The initial a z dynamics behavior (first 
plot) is characterized by a clear structure of beats, with beat 
and angular frequencies given by uii and f2 + uji/2. Due the 
decoherence processes, the modes related with the frequen- 
cies LOi evanesce in the characteristic time 1/T, destroying the 
beat structure (second plot), for which an oscillatory regime 
of same frequency Q, of the applied field arises. Here we as- 
sumed — 750. Solid and dashed curves represent the 
numeric and the analytic results obtained. 



tionary state, oscillating with frequency 2f2 (inset of fig- 
ure 21) and small amplitude around the asymptotic value 
a x (t — > oo). Here we observe that this asymptotic value 
changes dramatically when we leave the resonance con- 
dition. This feature can be understood as follows. At 
resonance, the driving field causes perfect transitions be- 
tween up and down eigenstates of a x , but because the 
system-bath coupling destroys the coherence along the 
process, we see that after a very long time, t ;» 1/r, 
a completely random emission-absorption process sets in 
and makes the expectation value of a x vanish. For the 
off-resonance case, the driving field no longer produces 
100% probability transitions in the system. Again, the 
coupling with the bath destroys the coherence of this 
process, and favours the occupation of the lowest energy 
state. Thus, as we are considering a small applied field, 
i.e, £o/A <C 1, the bath drives the system towards its 
thermodynamical equilibrium state at zero driving field, 
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Figure 7: The time evolution of the expectation value a x (t) 
using the approximate analytic solution Ea llll (dashed line) 
and the exact numeric calculation (solid line). Plot A presents 
the behavior for the resonance case f2 ~ A — 7. The oscillation 
frequency observed is u>i ft! eo. Plot B shows the case off 
resonance, f2 = 0.8A. The main oscillation frequency is u)%. 
The estimate relaxation time for these plots is 7\ = 12ps. 
The temperature used is T fa IK. 

a x — > tanh ( 2 k^T j ■ Since we have assumed HA IcbT, 

this implies tanh ( 2 ^ T ) ~ 1. How close the system will 

approach that asymptotic value depends on both the in- 
tensity and frequency of the driving field. 

All the features discussed so far were general conse- 
quences of the system-bath interaction, and do not de- 
pend on the specific form of the spectral density of the 
bath. Nevertheless, the first plots of figures [3] and [4] 
show a particular behavior for short times: two well de- 
fined regimes of oscillation frequencies can be observed. 
This phenomenon is a signature of the rich structure of 
the effective spectral density, Eq. ((3]). In fact, because 
J e ff(u) has poles in the complex plane, new modes of 
oscillation arise in the spin dynamics. Those new modes 
have lifetimes imposed by the imaginary part of the poles, 
whereas the natural frequencies of the system are shifted 
by their real part. As we could observe from J e ff(io) 
J3]), the most important modes ( those with long life- 
times and measurable amplitude of oscillations) due to 
the poles are associated with the peak of J e ff(u>), at 
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Figure 8: The time evolution of the expectation value a z (t) 
using the approximate analytic solution Eq[Tl] (dashed line) 
and the exact numeric calculation (solid line). Plot A presents 
the behavior for the resonance case fl ~ A — 7. uji and 
Q, — LUi/2 are the beat and the angular frequencies respec- 
tively. Plot B shows the case off resonance = 0.8A. We 
can identify approximately the beat and angular frequencies 
as loi and Q + lo\/2. The temperature used is T « IK. 



which a strong correlation between orbital and spin de- 
grees of freedom takes place. These modes, (see the 
first term of Eq. Ijlip ). have lifetimes determined by the 
width I of the bath resonance, and frequencies fl s — fiw, 
with = ft, = Q - ei and fl^ = ft + e x . 

The amplitude of oscillation, assuming a weak driving 
field, to /A <C 1, can be approximated by m^^y* (f° r 
A =^ Q s ). Therefore, if one experimentally reaches the 
conditions M^^k ~ 0(1) and £/A -C 1, it would be 
possible to determine some physical parameters of spin- 
bath coupling, e.g., the electron-phonon constant cou- 
pling S s . 

Figures [5] and [6] show the a z component for the same 
conditions previously discussed. As it can be noticed, 
there is a clear structure of beats in its dynamical be- 
havior. At resonance, figure the beat and angular fre- 
quencies are, respectively, given by cJi and fl + LUi/2. For 
this component we do not verify two oscillatory regimes 
for times of the order of l/£, This occurs because of the 
suppression of the pole's contribution by the strong os- 



cillations of the remaining terms. As a consequence of 
the decoherence process imposed by the bath of oscilla- 
tors, other important feature is that, for long times, the 
modes related with the frequencies u>i evanesce, destroy- 
ing the beat structure. These modes are associated with 
the coherent oscillations due to the driving field. Thus, 
after a long time, only an oscillatory regime of frequency 
Q remains. The characteristic time for that is of the or- 
der of 1/r. The regime slightly off resonance ft = A, 
figure [U practically does not change. In this case the 
beat and angular frequencies are given respectively by 
uji e Q — u>i/2. This shows that the component a z is not 
so sensitive to how close the external field frequency is 
to the natural resonance frequency of the spin. As it is 
directly coupled to the external field, this component has 
a stronger dependence on its amplitude. 

Now, we will focus on the specific case in which we 
use bulk physical parameters for typical quantum dots 
size. We assume here the piezoelectric electron-phonon 
interaction , s = 3, and quantum dots frequencies cl>o = 
15meV. The external fields applied to the dots are as- 
sumed to be such that ftA/fc^ = Q.7K and he^/kB = 
0.2K, and the Debye frequency hu>u/kB — 360i^. Fig- 
ures [7] and [5] present the a x and a z dynamics using the 
bulk physical parameters 13 = 3000m/s, m* — 0.063m e 
and 63 = 355. The temperature used for the calculation 
is T « IK. 

The plots A in figures [7] and [8] show the dynamics 
of the expectation values a x (t) and a z (t), respectively, 
for the case close to the resonance, fl ~ A — 7. We 
can see that the dynamics of a x has a well defined fre- 
quency of oscillation, given by u>i ~ eo, with amplitude 
close to |<Ta;(0)|. Because of the weak coupling limit, the 
terms in Eq. fTTj) due to the poles of the spectral density 
are negligible compared to the others. Therefore, this 
explains why we do not distinguish between several dif- 
ferent regimes of relaxation. The relaxation time T\ is 
estimated to be of the order of 1/r = 12/iS. For the <j z 
component we see the expected beat structure, with well 
characterized beat and angular frequencies given by uji 
and £1 — Wi/2, respectively. 

The plots B in figures [7] and [8] present the dynamics 
for a case off resonance, fl = 0.8A. Now we can see that 
several frequencies contribute to the spin dynamics. For 
a x , the main oscillation frequency is still w\, however the 
amplitude of oscillations has decreased more than one 
order of magnitude. The a z dynamics has also changed. 
There is still a structure of beats but more frequencies 
appear in the spectral decomposition. uj\ and fi + wi/2 
are the beat and angular frequencies for this case. For 
the very long time regime, t ^> 1/r, the component 
a x reaches the stationary equilibrium oscillating with 
frequency 20, and very small amplitude G, around the 
value a x (oo). For a z we observe that the beat structure 
disappears, giving place to oscillations of frequency Q 
and constant amplitude ^ f (eo, Q)(a x ((X)) + G/2). 

In conclusion we have obtained approximate analytic 
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solutions for the Blcoh-Redfield equations of dissipative 
spins in a driving field. The complex structure presented 
by the effective bath spectral density brings a richer 
dynamical behavior for the spin trapped in a quantum 
dot. We showed that the existence of the peak in the 
spectral density introduces a new temporal scale in the 
problem, l/£, the peak's width, which clearly separates 
distinct time evolutions of the component a x . The 
first one, short times, t <C l/£, is basically determined 
by the characteristics of the peak. We saw that the 
oscillations presented in this regime have their natural 
frequencies shifted by the bath resonance, Sl s , and their 
lifetimes are determined by l/l. Thus, in principle, 
this regime could be used to infer intrinsic quantities of 
the system-bath coupling , such as the electron-phonon 
constant S s . On the other hand, the long time regime, 
t ^> l/£, was shown to be governed by the external 
fields. We found that the relaxation rate T for this 
regime presents a strong peak for frequencies close to the 
bath resonance, revealing the signature of the effective 
spectral density. In addition, we were able to identify the 
shift induced by the system-bath coupling in the natural 
frequencies of the system. We found that this shift, 



7, has the form of a temperature dependent Lamb shift^ 
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VI. APPENDIX 

In this section we present the explicit form for the 
Bloch-Redfield coefficients used in the main text. 

The first step is to construct the matrix elements of 
the non-dissipative time evolution operator U(t,t'), 



lmU RR (t, t')U RL (t, = ( sin[(f - t')Sl] + sin[(t + t')Sl](l - cos[(t - «>i])) 



(1 + sm(d>)) 2 . ,. .. .„ X1 (1 — sin(0)) 2 . ,, 

"- sm[(t - t')(Q + ei)] + i w " sm[(i - t')(Sl - ej)], 



ReU RR (t,t')U RL (t,t') = 
\U RR (t,t')\ 2 -\U RL (t,t')\ 2 = 



8 

^^sin(t'0)sin[(t-i')ei] 

COS 2 (0) 



cos(t'n)(l - cos[(t - t')ei]), 



sin(20) 



cos[(t - t')Q] + cos[(i + t')Q}(l - cos[(< - t')e 1 }) 
y " cos[(i - If) (SI + ei)] - ^ — — — ^— cos[(t - t')(Sl - ex)}. 



(15) 



r 



If the limit ^p-\uj^o exists, and J (to) does not 
have poles in the real axis, we can write the 
real and imaginary parts of M(t) as Re[M(t)] = 
i r o °° du J(uj) cos(ujt) coth(/3fiw/2) and lm[M(t)] = 
— i J °° duiJ(uj) s'm(ujt), respectively. Now we are in po- 
sition to calculate the coefficients Ai(t) and Tij(t), As 



pointed out previously, if we first perform the time in- 
tegrals, we can reach very useful expressions for the fre- 
quency analysis of our problem. Following this proce- 
dure, we obtain the fluctuation coefficients A x and A y 
and the rate T xx (t) as (with analogous expressions for 
T xz (t), and T yz (t)) 



A x (t) = - 



l-cos(2fit) )l[ 0) [Sl] J(x),t] - 



cos 2 (</>) 

cos 2 (0) cos(2fit) + (1 - sm((/>)) 2 

47T 



cos 2 (0)cos(2m) + (l + sin(0)) 2 ^ (0) 
— I 1\ + er, J(x),t 



l[ 0) [Sl-ex;J(x),t] 



C -^ si n(2Slt)(l^[Sl-,J( X ),t] 



2n 



- \ (4 0) [SI + ex ; J(x) , t] + 4 0) [SI - ex ; J(x) , t] ) 



(16) 
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A y (t) 



sin(20) / 
cos(0) 



2tt 
cos(0) 

C0S 2 (</)) 



sm(nt)l[ Q) [ft; J(x),t] + cos(m)4 0) [f2; J(x),i 

sin(Ot) (\l + sin^))/} 03 [fi + ei; J(x), t] - (1 - sin(0))/j O) [O - e i; J(x), t] 
cos(n<)( (1 + sin(0))4 O) [fi + ei; J(a;),<] - (1 - sin(<£))4 0) [fl - er, J(x), t] 



2?r 



l + cos(2Qi) ]J 3 [fi,T,t] -sin(2fii)( / 4 [fi,T, t] - -(J 4 [fi + ei,T,t] + - ei,T,t}) 



(17) 
(18) 



cos 2 (0) cos(20<) - (1 + sin(c/>)) 2 



4tt 

where we have defined the integrals 

h[y,t] = i 

h[y,t] = 

h[y,T,t] = 

h[y,T,t] = 



7 3 [n + ci,T,t]- 



cos 2 (</>) cos(2Qi) - (1 - sin^)) 2 



4-7T 



/ 3 [n-ei,T,t], 



dujJ(ui) 
dujJ(uj) 



sin((a> — y)t) sm((u> + y)t) 
ui-y uj + y 

sin 2 ((w -yU) sin 2 ((u + y)f) 



dujJ(ui) coth 
dujJ(w) coth 



w -y uj + y J 

( Hlj \ /sin((w — y)t) sm((uj + y)t) 



(19) 



\2k B T 

fku 
2k B T 



u>-y 



uj + y 



sin 2 ((w-y)|) sin 2 (( W + y)|) 



uj-y 



uj + y 



(20) 



I 

The special form of the integrals /„ allows us to find fies the conditions J (0) = J (oo) — > 0, we can write for 
approximate analytic expressions, which are asymptot- t 3> l/i~2 c , 
ically exact for i — > oo. Considering that J (oj) satis- 



A(y,i)«|j(y) + ^R<'s 



J(«o) 



,it(z —y) 



TV 



h(y,t) wi jf C 4;J(w)P 



1 



Res 



1 



o it{z +yY 



zo + y 



- Res 
Res 



-it(zo-y) 



zo-y 

e —it(z +y) 



u>-y u+y 



ITT . 

-y^Res 



J(2o) 



o 2it(z -t/) 



/7T 



Res 



J(zo) 



zo-y 

,2it(z +y) 

zo + y 



za + y 

- Res 

— Res 



J(z ) 
J(z ) 



-2it(z -y) 



zo-y 

e -2it(z +y) 



zo + y 



(21) 



(22) 



and similar expressions for I 3 and h- The residues in the 
above expressions are calculated for the poles of J (cj) 
inside of the semi-circle centered at z — fl c and radius 
fl,r^. In this calculation, we have assumed that y is not 
a pole of the function J (uj) . Here we have replaced the 
upper limit of integration by the cutoff frequency fl c - 

Expressions f2Tj) and (|22| show us that the time- 
dependent terms of /„ have lifetimes determined by the 
imaginary part of the poles of J (uj) in the complex plane. 



In our case, the poles of J (uj) correspond to those of 
J e ff (uj) defined by EqEl 

Figure [9] sketches the I\ and I2 behaviors assuming 
J(u>) as an Ohmic function J(u>) — O.OSwe -20 " and the 
spectral density J e ff(u>) (Eq. [3]). Comparing the results 
using the analytic form Eqs. [2l1f22l and the exact numeric 
calculation, we find a good agreement for t 3> 1/Vd, 
specially for the case g(u>) = J e //(w) where the terms 
due the poles dominate the short time regime, t <l/£. 
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